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Abstract 

We discuss persistent currents for particles with internal degrees of freedom. 
The currents arise because of winding properties essential for the chaotic 
motion of the particles in a confined geometry. The currents do not change 
the particle concentrations or thermodynamics, similar to the skipping orbits 
in a magnetic field. 
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1. Introduction 

Windings of random walk trajectories are of considerable interest in vari- 
ous fields of condensed matter physics, from polymers and DNA to supercon- 
ductors and Bose- Einstein condensates (see e.g. and references therein). 
In this paper, we consider the influence of winding on the internal degrees of 
freedom of particles dynamics. 

Note that winding dependent degrees of freedom are not exotic. The 
most obvious of these is, of course, the spin, which is of particular interest 
in view of the recent spintronic surge of activity. Because of the spin-orbit 
interaction, the spin is directly related to the winding of particle trajectories. 
Accordingly, our analysis is intimately related to the so-called spin Hall effect 
in electron transport as well as to other types of spin transport, like, 



for example, photon (light) propagation in disordered media (since photons 
also have non-zero spin). 

Other degrees of freedom, both classical and quantum, can also be wind- 
ing dependent, as e.g. geometric phases (the Berry phase). Thus, the results 
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of the present work are applicable not only to spin transport but also to any 
situation that can be reduced to random trajectories in a confined geometry 
with the dependence of certain characteristics on winding properties. 

The paper is organized as follows. Section 2 presents the concept of wind- 
ing and the corresponding terminology. Section 3 deals with the spin degree 
of freedom as directly dependent on a winding angle. Section 4 contains nu- 
merical simulations in the case of diffusion in a rectangular 2d domain with 
periodic boundary conditions in the longitudinal direction and hard walls 
in the transverse direction. Section 5 proposes a simple theoretical model 
of classical random walks with so-called "soft walls", where the geometric 
constraints are described by an external confining potential in the transverse 
direction. Section 6 contains a brief discussion of the results. 



2. Winding of curved trajectories 

Let us introduce the winding angle 6 as the vector sum of all the turnings 
of a particle trajectory x(X)- The angle is called the total curvature in the 
literature (see e.g. 11, lUl). It is given by the integral of the differential 
curvature over a natural parameter (the local time of the trajectory). The 
differential curvature n of the trajectory is defined as 

where the symbol " ' " denotes differentiation with respect to t and the 



symbol x denotes the vector product. Correspondingly, 6 has the form [13 



/"* f'^ V X V 

0=1 hi{t)\x{t)\dt = j —dt, v = x. (1) 
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Note that the total curvature of a closed planar trajectory, equal to the 
"winding number" or "turning number" multiplied by 27r, is a homotopy 
invariant. 

It is worth mentioning that the term "winding angle" is also used in 
polymer physics (see e.g. |3|) where, however, it rather refers to 

/■* a; X i , 



Thus, is the winding angle of the radius vector of the trajectory in config- 
uration space, whereas yields an analogous quantity but in velocity space. 
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The distribution of i^'s for Brownian trajectories has been studied for a 
long time (see e.g. j2|, [s], [13]), while that of 0's has not been analyzed so far. 
On the other hand, it seems that is a more appropriate quantity for certain 
problems, for example those including spin (spin Hall effect, etc), as we will 
argue below. 

We consider for simplicity the 2d case. Generalizations to 3d and tensors 
are fairly simple, but the resulting formulas are more involved and less intu- 
itive. In the planar case the winding angle vector is normal to the plane, 
hence it suffices to consider 9, the unique non-zero component of d. 

The transport of the winding angle is described by the average current 
density 

3{x) = {ev5{x-x{t))), (2) 

where the symbol (...) denotes averaging over a set of Brownian trajectories, 
or, in the discrete case, of random walks. 

In the case of an infinite plane (or a plane with periodic boundary con- 
ditions) and of equilibrium, the particle density does not depend on the 
coordinates and the average (6*) vanishes identically as well as the average 
of the winding current density If, however, a confined planar geometry 
is considered, then {6) will still be vanishing but, as we will show below, 
this will not be the case for the winding current, because of the geometrical 
constraints. Namely, due to the mere existence of boundaries, there will exist 
trajectories whose contribution to the winding will not be canceled by that of 
other trajectories. As a result, there will necessarily appear an anti-parallel 
edge current with ^ > and ^ < and j will not identically vanish anymore. 
These surviving edge currents can be viewed as analogs of spin edge currents 



in the spin Hall effect [15 



3. Spin-orbit interaction and winding 

Before considering the winding currents, it is useful to discuss briefly an 
important example of spin-orbit interaction where the spin degree of freedom 
depends on the winding angle 9. 

It seems known that in electronic transport in the presence of charged 
impurities the electron spin depends on the winding of their trajectories. 
This dependence manifests itself, for instance, as the so called skew scattering 
in the extrinsic" spin Hall effect 0, [sj. However, in the literature on the 
subject, the concept of trajectory winding is mostly implicit. 
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Recall that in electronic transport impurities interact with electrons and 
thus the scattering angle, hence the deviation of the electron trajectory from 
a ballistic straight line, depends, in relativistic theory, on the spin. More 
generally, it is known that if, for any reason (not necessarily related to elec- 
trons scattering on charged impurities) the trajectory of a relativistic particle 
is not a straight line, the curving of the trajectory changes the spin degree 
of freedom of the particle. This purely kinematic relativistic effect is due to 
the successive non-parallel Lorentz transformations (boosts) that lead to the 
so-called Wigner rotation of the spin, or the Thomas- Wigner precession. 

The variation of the precession angle of a classical rotating moment re- 
sulting from an infinitesimal turning of the velocity has been obtained by 



Thomas 16-18 



dip ^ -^v X dv , (3) 

where is a vector of the precession angle, v is the velocity, dv its variation, 
and c is the speed of light. 

Eq- © gives the leading relativistic approximation for the precession 
angle. It is worth mentioning that the exact relativistic expression is still 
a matter of controversy (see the recent articles 19-2l|). Note also that the 
Dirac equation allows for a more precise study of this effect. However the 
classical result (E]) of Thomas is sufficiently accurate to describe well the spin- 



orbit splitting in atoms [22|. It has played an important role in the genesis 
of the quantum theory of radiation (see e.g 18|). 

Comparing ([1]) and (|3]), we find a simple relation between the moment 
precession angle and the winding angle of the trajectory 

dip ^ ^de , (4) 

where dO is the trajectory turning angle. Note that, in many cases, the pre- 
cession angle during the observation time is essentially equal to the winding 
angle multiplied by a small parameter v'^/c^ (e.g. the case of fermions at 
temperatures below Fermi temperature, or any particles in the case of elastic 
scattering) . 

Let us now argue that the spatial curvature of a trajectory is equivalent 
to the action of a magnetic field on a spin degree of freedom. Indeed, it is 
known that an homogeneous magnetic field in the rest frame of a classical 



spin results in its Larmor precession [23 



sx H 

2mc 
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where s is the spin in angular momentum units, H is the magnetic field, g 
is the gyro-magnetic factor, and m is the mass. Thus, an effective magnetic 
field with a Larmor frequency equal to the frequency ip of Thomas- Wigner 
precession {s = s x cp) has the form 

2mc m . 2mf ^ • . . 

ileff = <^ ~ V X V = 6 . (5) 

ge egc egc 

The spin-field interaction energy is then 

2mc 2c2 c2 

Hence, the direction and the magnitude of the effective magnetic field are 
intimately related to the winding angle. 

Returning to random walks of electrons in soli ds ( which can result from 
electrons scattering on charged impurities, see e.g. |2J|), we note that in this 
case it follows from that the effective magnetic field of the Lorentz trans- 
formation of the impurity electric field can be included in the Wigner-Thomas 
precession and vice versa, thus forming a complete spin-orbit interaction. 

We also note that, leaving aside electrons in solids, the above discussion 
can refer in general to the diffusion of any particles with non-zero spin or 
magnetic moment in disordered media, for example spin 1 photons whose 
motion in a medium with random refractive index can be described as a 
diffusion process 24j-|26|). as well as in foams, metamaterials, dense plasma, 
etc... (see e.g. |27| and references therein). 



4. Numerical simulations 

Let us now turn to numerical simulations which do support the existence 
of winding edge currents. We consider the strip of width Li along the hori- 
zontal axis and of length L2 along the vertical axis on a planar square lattice 
of period 1. In the longitudinal direction periodic boundary conditions are 
imposed at X2 = 0, L2. In the transverse direction reflecting (hard) walls 
enclose the strip in a flnite width xi = 0, Li. 

We consider Np independent particles whose trajectories are random 
walks on the lattice. The diffusion time step r is assumed to be unity and 
random hoppings are allowed only to the nearest neighbor sites. Let 

k 

x{k,i) = Ax{m, i) 

m=0 
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and 6{k,i) be respectively the position and the winding angle of the ith 
particle after k jumps. It is clear that each next turn (left or right) increases 
or decreases 9{k,i) by n/2. We assume, for the sake of definiteness, that a 
backward scattering does not change the winding angle of the trajectorjQ. 
The time averaged winding current j{x) at a point x of the lattice is after 
steps (i.e. during the observation time t = N) 

1 ^ 

= ]V 5^ 5Z > (6) 

i k=0 

where Sj-y is the Kronecker delta. 

We analyzed a wide range of parameters Li, L2, N, and Np. Typical 
results are shown in FigH] for Np = 4000, = 500000, and strip sizes 
Li = L2 = 100. To display the 2d vector field j, we use a 2d palette as in 



Fig, 1(b) The palette center corresponds to zero vector field and its vertical 
and horizontal edges correspond to the range of the components ji and j2 of 
j{x), respectively. The palette is divided in four sectors of different colors 
with a radial variation of the color intensity. This allows to distinguish the 
approximate direction and magnitude of the current. 



The numerical simulations for j{x) are shown in Fig 1(a) In the bulk of 
the sample only small-scale fluctuations of the winding current can be seen. 
However, things differ strongly from the bulk in the narrow neighborhoods 
of the walls: winding currents do exist along the edges. They materialize as 
almost continuous and narrow lines on the left and right edges of the sample 
(blue and red, according to the direction). The non vanishing longitudinal 



component of the winding current, averaged over X2, is displayed in Fig 1(c) 
It is concentrated within one lattice site from the edges and directed down- 
ward near the left edge and upward near the right edge. 

Fig 12] displays the absolute value of the current density \j{x)\ and how it 
concentrates on the edges with an increase of the simulation time A^. 

Since the system is in equilibrium, i.e., external fields are absent and the 
particles are distributed uniformly, winding currents can only originate from 
the confining boundaries and the chaotic kinetic energy of the random walk 
trajectories. Thus, these currents are by nature persistent. 



^We did consider other variants, but they do not modify qualitatively the effect, though 
it is rather obvious. 
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5. A simple solvable model 



Let us consider a simple model where tlie effect of winding currents is 
explicit. It consists of planar continuous Brownian trajectories described by 
an Ornstein-Ulilenbeck process in a potential U applied along the xi axis 
and infinitely growing as |xi| — )■ oo. The potential plays the role of the "so/^ 
walls" confining the Brownian particles, and thereby encodes the geometric 
constraints. As in the previous section, periodic boundary conditions at 
X2 = 0, L2 are assumed in the longitudinal direction. 

The stochastic (Langevin) equations of motion are 

Xi = Vi (7) 

Vi = 7^ PVi + Q 

m oxi 

where x = (xi, X2) is the coordinate of the particle, v = (f 1, f 2) is its velocity, 
m is the mass of the particle, f3 is the viscosity and ^ = (^1, ^2) is the random 
force assumed to be the Gaussian white noise, i.e. 

m)) = 0, m%in) = - n- (s) 

Since we are interested in the winding of trajectories of the above random 
process, we consider the joint probability density of the coordinate, velocity 
and winding angle 9 

pit, X, V, 6) = {5{x - x{t))S{v - v{t))S{e - e{t))) . (9) 

The 6 integration of ([9]) yields the density P(t, x, v) satisfying the usual 
Kramers equation jisi (Fokker-Planck equation with external potential) 

9 ^ 1 dU d ^ ^ d „ 
-P = -v,—P + - ——P + P—ViP+--^P. 10 

ot oxi m oxi ovi ovi 2 ovf 

where D is defined in and the standard summation convention over re- 
peated indices is understood. 

It is easy to check that the probability density P in ffTOj) converges as 
t — 00 to a stationary Maxwell-Boltzmann distribution in xi and v with 
a normalization factor in the longitudinal direction X2 (see f|T5l) ). Note 
that one could also use reflecting boundary conditions or a weak confining 
potential in this direction. 
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The Fokker-Planck equation for p{t, x, v, 9) can be obtained as follows. 
Let us denote 

X> = 5{x- x{t))5{v - v{t))5{e - e{t)), 

then p = {Cl) and 



d_ 
di 



Then eq. implies 



d 



d 



d n dU 



dU d 



mv'^ ^ dxj dd^ 



d 



d 



m- (11) 



To deal with the correlation functions (^jO), we use the Furutsu-Novikov 
formula (see e.g. 29|, Chapter 10 and [30], Chapter 5) valid for any functional 
Rt[^] of the 5-correlated Gaussian random functions .^i and .^2 



2\50(t) 

We obtain after integration by parts 

= ^ 9 /6x,it)^^ _ D^/6v^^^^ _ D d / 66(1) 



2 dxk \ 6^,{t) 



2 dvk\S^j{t)^/ 2 de\6^j{t)^ 



:i2) 



To find the variational derivatives of Xk and Vk with respect to ^j, we apply 
the operation 6/6C,{t) to the integrated Langevin equations, then use the 
principle of causality and take the limit t ^ t. This yields 



t 



> 



and 
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The variational derivative of 9 in f ll2l) is calculated using ([T]) 



— eik I dti- 





and finally 

t 

66{t) _^ ^ Vi{T) eik /"^^ r T "^^^^ _e 



v^{t) m J 



T 

Thus, the correlation functions (^jf ) are expressed via the partial derivatives 
of p 

,^ , , , D d D Vi d 

We assume for the sake simplicity that the diffusion coefficient D does not 
depend on the velocity and we take into account the relation eijejk = —Sik- 
We obtain the Fokker-Plank equation 

dt^ dxi^ dvi\mdxi 2 dvf^ 

l_ v^dJ^d_ Vj d'^ D 

In the case at hand where the potential U depends only on the transverse 
coordinate Xi (modeling the soft walls of a deep valley) the Fokker-Planck 
equation becomes 



d d d 1 d D d'^ 

T^P + Vi^P--^ U'{xi)p-l3—Vip-——^p (13) 

at oxi ovi m ovi 2 ai; 



m f ^ 06 OV206 f ^ ovioa 2 06^ 

Eq. f fT3|) differs from the Kramers equation by the terms on the right-hand 
side. The physical meanings of the first and fourth terms are the most 
interesting. It is already clear here and will be even more clear below that 
the former is the current in the X2 direction and the latter is responsible for 
the diffusive spreading of 6. 
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We also note that after integration over 9 fll3p coincides with the Kramers 
equation (fTUil . 

Eq. f|T3l) allows us to study the transfer of arbitrary 6'-dependent quan- 
tities as well as non-equilibrium situations corresponding to arbitrary initial 
conditions. In the particular case of the winding current ([2]), which is lin- 
ear in 6', the steady-state current can be found without actually solving the 
equation. 

Indeed, it is clear that the transverse component ji of the mean current 
density j vanishes. To find the average longitudinal component j2 let us 
multiply (fT3|) by Ov2 and integrate over v and 6. We obtain after integration 
by parts 

^ = -/3j2 + -U'{x,) [ dv%P{t, X, v) . (14) 
ot m J 

The probability density P{t, x, v) approaches for large time the stationary 
distribution 

P I [i ^ 2[i 



"•"-^rA-!^'-'-^''^'''')' '''' 

where 

is a normalization constant. It follows then from fll4p and fllSp that for 
t» the current 22 relaxes to its stationary form 

J2(xi) = i^^^\^^)Pb{x,) (16) 

where 

is the steady state (Boltzmann) probability distribution in the transverse 
direction. 

We obtain finally for the local average density of winding current per unit 
length in the longitudinal direction 

where n = Np/ L2 is the number of particles per unit length in the longitudinal 
direction. 
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The average winding current density ( II 7p is proportional to the gradient 
of the potential and to the local density of diffusing particles, which is deter- 
mined by the temperature and the height of the potential at a given point 
Xi. One can easily see that the currents are antiparallel on opposite walls. 
The gradient of the potential and the Boltzmann distribution also determine 
the currents width. Particular examples of currents profiles are shown in 
Figure E] for an harmonic potential U{xi) = x\ and a trapezoidal potential 
U{xi) = ?7(— 1 — + 1| + ri{xi — l)\xi — 1|, where ri{x) is the Heaviside 

step function [ ri{x) = 1 for x > 0, and r]{x) = for x < ]. 

Note that formula (fT6l) can be interpreted (up to a small relativistic fac- 
tor) as describing the diffusion of magnetic moments in an effective magnetic 
field H (X [v X VU], thereby corresponding to the setting of the spin Hall 
effect in the potential U, as already mentioned above. We recall that in our 
approach there is no interaction, and no magnetic (or any other type) mo- 
ments. We are only dealing with distributions of trajectories winding, hence 
the above effect can be viewed as purely dynamical ("winding Hall effect"). 

Assume finally without loss of generality that U (0) = 0. Then the spatial 
average of J2 in flTTl) over (0, 00) (or (—00, 0)) is 



This can be viewed as a manifestation of the weak dependence of our results 
on the concrete form of the confining potential (considering that n/Z is the 
effective planar density of particles). 

6. Conclusion 

We demonstrated that random walks or continuous Brownian motions 
in a confined planar geometry generate dissipationless persistent currents of 
degrees of freedom associated with the winding of their trajectories. This is 
a purely edge effect where the currents are formed in a close neighborhood of 
the boundaries. In the case of the "soft walls" , where the geometric confine- 
ment is provided by an infinitely growing transverse potential, the channel 
width of the current is determined by the gradient of the potential and the 
thermal energy of the particles. In the case of reflecting boundaries ( "hard- 
walls"), the channel width is just a single diffusion step, according to the 
numerical simulations. 



/ 

Jo 



00 




nD 



(18) 
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In the microscopic world, the spin and geometric phase are examples of 
degrees of freedom which depend on the winding. Note that it is commonly 
believed that in the spin Hall effect the dissipationless edge currents arise 
in the conditions of the quantum 'Hntrinsic" regime, where dissipation and 
chaos are absent because of quantization. Since, however, the above "eodrin- 
sic" spin Hall effect corresponds to an ordinary diffusion, thus to a sufficiently 
high dissipation, one would think that persistent spin currents should not oc- 
cur. However, as we have shown, this is not the case simply because of the 
chaotic kinetic energy of the random trajectories of the particles. 

Similar effects may occur in purely mechanical systems where relativistic 
or quantum phenomena are absent. For example, the torque of a compass 
arrow in a moving car is also a degree of freedom which " feels" the curvature 
of the car trajectory. Thus, the effect studied here can be of interest beyond 
the scope of phenomena associated with the spin Hall effect. 

In the case of the Ornstein-Uhlenbeck model with soft walls we have 
obtained an analytical expression for the persistent current of the winding 
via the Boltzmann distribution in an arbitrary transverse confining potential. 
The kinetic equation takes into account the winding of particles and, in 
principle, allows to explore the joint probability distributions in the single- 
particle phase space of the degrees of freedom with an arbitrary dependence 
on the winding. 
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LIST OF FIFGURES 



Figure 1. (a) - the winding current j. (b) - the 2d palette, (c) - the 
winding current component j2 averaged over X2 as a function of Xi. 

Figure 2. The modulus for different values of with Np — 3000. 

The darker the dot, the bigger the modulus of j. 

Figure 3. The winding current density in the harmonic and trapezoidal 
potentials and the profile J2{xi) in arbitrary dimensionless units. On the 
profile charts, the left y-axis is J2 and the right y-axis is U. 
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Figure 1: (a) - the winding current j. |(b)| - the 2d palette, (c) ~ the winding current 
component j2 averaged over X2 as a function oi xi. 
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Figure 3: The winding current density in the harmonic and trapezoidal potentials and the 
profile J2{xi) in arbitrary dimensionless units. On the profile charts, the left y-axis is J2 
and the right y-axis is U . 
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